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The tight binding model for an electron on an anisotropic 
triangular lattice in a uniform magnetic field is studied using 
a decimation scheme. The model exhibits a transition from 
critical to localized phase and the phase diagram is described 
in terms of three nontrivial renormalization fixed points for 
the band edges. These subcritical, critical, and supercritical 
universality classes also describe the corresponding states of 
the quantum Ising chain in a modulating transverse field. The 
only exception is the conformally invariant point of the Ising 
model which has no analog in the electron problem. 
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The problem of a two-dimensional Bloch electron in 
an irrational magnetic flux is important both theoreti- 
cally and experimentally. In this paper, we study the 
electron on a triangular lattice with anisotropic couplings 
between the nearest neighbor sites. We use a renormal- 
ization approach to obtain the phase diagram and uni- 
versal scaling properties of the model. The results are 
compared with those of the Ising model in a transverse 
modulating field of periodicity incommensurate with the 
periodicity of the spin chain. With the exception of con- 
formally invariant point of the Ising model Q, the 
band edge as well as the band center states of the two 
models belong to the same universality class. 

In the tight binding approximation, the energy spec- 
trum for a single Bloch band of an anisotropic triangular 
lattice with couplings t\ and t% is |Q| , 

e(k x ,k y ) = -t lC os[^V3k x -k y )] 

— t\ cos[ )] + t% cos(ak y ) (1) 

where a is the lattice spacing. Using Peierl's substitution 
hk — > p — -A with p = jV and the Landau gauge A = 
(0, x, 0)H for the vector potential A, the two-dimensional 
problem is reduced to the one-dimensional tight binding 
model (TBM) f§ 

cos[vr(icr + (<f> + ^r))]tpi+i + cos[7r(i<7 + (0 - 



+ Acos[27r(icr + 4>)]ij)i = Eipi 



(2) 



which we call the triangular model. Here, i = 2x/(ay/3) 
and a is the dimensionless magnetic flux in units of Bohr 
magneton, a = (eH/hc)a 2 and A = ti/ta- The isotropic 
triangular lattice case t\ = ti was studied by Claro and 
Wannier [Q where the model was found to exhibit the 
hierarchial structure in the energy spectrum similar to 
that of the square lattice case studied by Hofstadter H . 

In this paper, we show that as a function of A, the 
model exhibits a phase transition from critical (C) to lo- 
calized (L) states. This is in contrast with the case of 
the square lattice where the anisotropic model exhibits a 
transition from extended (E) to localized states. There- 
fore, unlike the square lattice case, the model does not ex- 
hibit the E phase of KAM type and hence the weak cou- 
pling limit of the model is nontrivial. Furthermore, the 
critical phase exists in a finite parameter interval. We use 
our recently developed decimation method ||,|6| to study 
the scaling properties in various parameter ranges. At 
the band edges, the universality classes of the model are 
described by the following three non-trival fixed points 
of the renormalization operator (the term "fixed point" 
is used also for periodic cycles of the renormalization op- 
erator as the latter are fixed points of an higher-order 
iteration of the operator): 1) the weak coupling (sub- 
critical) fixed point (A < 1), 2) the critical fixed point 
corresponding to A = 1 describing the onset of the C- 
L transition, and 3) the strong coupling (supercritical) 
fixed point (A > 1) describing the fluctuations in the ex- 
ponential decay in a localized wave function JtJ . 

The C-L transition and the universal properties of the 
triangular model (2) are compared with those of the 
quantum Ising model (QIM) in a transverse field hi for 
which the Hamiltonian reads 



ff = -£[<r?Vf +1 + 2/M7f 



(3) 



Using the methods described by Lieb et al. |§ , the eigen- 
value equation for the spin problem can be written in a 
TBM form, which for hi = Acos[7r(icr + </>)] becomes 
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cos[7r(z(T + (j>)]ipi+i + cos[7r(icr - 
+ Acos[27r(icr + (f))]ipi — Eipi 



4>M 



(4) 



Here E = -A + (E 2 /4 - 1)/(2A), where E is the energy 
of the Ising model (3). Q 

The QIM was recently studied for a equal to the in- 
verse golden mean a = (75- l)/2. The model exhibits 
critical states for A < 1 where A = 1 is the localization 
transition point. In the QIM, this transition has an addi- 
tional significance as it also corresponds to the magnetic 
transition to the long range order. At this point, the low- 
est energy quantum state is believed to be conformally 
invariant. [g| Our recent decimation studies || confirmed 
this phase diagram of the model and showed that the sub- 
critical regime was described by a unique fixed point of 
the renormalization flow. However, at the onset of lo- 
calization, the renormalization flow was attracted by a 
different fixed point resulting in a different universality 
class. 

It is rather interesting to note the similarities between 
these two seemingly unrelated problems. Both are de- 
scribed by a nearest-neighbor TBM where the diagonal 
as well as the off-diagonal terms are modulating: the pe- 
riodicity of the diagonal term is twice the periodicity of 
the off-diagonal term. The only difference between these 
two models is the fact that unlike the QIM, the triangu- 
lar model shows a relative phase difference between the 
diagonal and off-diagonal term which is a dependent. 

We applied the renormalization approach || to the tri- 
angular model. Our main focuss is to compare the phase 
diagram and the universal properties of the triangular 
model with those of the QIM. The TBM is written in the 
following form, where all sites of the lattice except those 
labelled by the Fibonacci numbers F n are decimated: 



fn(i)i>(* + F n+ i) = + F n ) + e n (i)i)(i). 



(5) 



The additive property F n +i = F n + F n -i provides exact 
recursion relations for the decimation functions e„ and 

Ae n (i) 



1 + Af n (i) 



(6) 



, M fn-i(i + F n )f n (i + F n ) 

A = e n _i(z + F n ) + / n _ x (i + F n )e n (i + F n ). 

In this renormalization approach, the C phase manifests 
itself in a nontrivial limit cycle of period p, which is typi- 
cally a multiple of 3, for the absolute values of decimation 
functions e n and /„ as n — ► oo. This results in nontrivial 
scaling ratios 



Ci = Hm \ip(F i+pn )/iP(0)\. 



(8) 



Here it is assumed that the phase <f> is chosen so that the 
main peak of the wave function lies at i = 0. In contrast, 
the scaling ratios approach unity in the E phase. 



In the L phase, we write the wave function as 

A = e-^ m (9) 

where 7 is the Lyapunov exponent which vanishes in the 
E and C phase, rj describes the fluctuations in the expo- 
nentially decaying part of the wave function. The spatial 
dependence of 77 is given by a TBM |7j which is studied 
by the decimation scheme described above. For the tri- 
angular lattice, 7 = ln(X) Q. Due to lack of duality, an 
analogous expression has not been derived analytically 
for the QIM. However, a numerical computation of the 
Lyapunov exponent [[lOj is found to be in agreement with 
this formula. 

Detailed study of the renormalization flow shows that 
for the triangular model (2), the decimation functions 
approach asymptotically the same 3-cycle for both the 
upper and the lower band edge (maximum and minimum 
energy). For the QIM, this identity between the scaling 
properties of the upper and lower band edge is true only 
in the subcritical and supercritical regimes but not at 
the critical point (see next paragraph). Furthermore, in 
both subcritical and supercritical regime, the triangular 
model and the QIM are described by the same univer- 
sality class, i.e. the scaling ratios £ for both models are 
identical. In the triangular model the wave function has 
exact symmetry at the band edges whereas in the QIM 
the wave function is only asymptotically symmetric. 

At the critical point A = 1 corresponding to the onset 
of localization, the universal characteristics of the trian- 
gular lattice band edges are the same as those of the 
upper band edge of the QIM. However, the lower band 
edge of the QIM has zero energy and the corresponding 
state is believed to be conformally invariant. Unlike the 
upper band edge 3-cycle, the renormalization flow for the 
decimation functions at E = converges to a period- 1 
fixed point. In addition, the wave function at the confor- 
mal point is asymmetrical vanishing on one side of main 
peak. H The " conformal" renormalization fixed point of 
the QIM does not map to any quantum state of the tri- 
angular model. 

Fig. 1 shows the self-similar wave functions in the three 
universality classes for the triangular model at the up- 
per band edge. They also describe the asymptotic wave 
functions at the upper band edge of the QIM. The dom- 
inant peaks in the wave functions are labelled by the 
rational approximants of a 3 and its harmonics. || The 
height of these peaks with respect to the central peak de- 
fine universal scaling ratios. For example, the universal 
scaling ratios of the peaks at three successive Fibonacci 
numbers in the subcritical, critical, and supercritical 
phases are (0.825,0.908,0.836), (0.238,0.303,0.291), and 
(0.267,0.311,0.121), respectively (see Fig. 1). These 
scaling ratios are identical to the corresponding scaling 
ratios at the lower band edge with the exception of the 
the conformal point in the QIM characterized by the fixed 
point scaling 0.415. 
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The study of the band center indicates how the rela- 
tive phase difference in the off-diagonal terms of the tri- 
angular model and the Q1M can lead to subtle symmetry 
properties in the asymptotic renormalization dynamics 
of the two models. At the band center with A = f , the 
wave functions are asymmetric in both systems leading 
to different 6-cycles and scaling ratios (also different from 
those characterizing the band edges) on the positive and 
negative side of the lattice. The scaling properties on 
the positive (negative) side of the triangular model are 
asymptotically the same as those of the QIM on the nega- 
tive (positive) side. However, the correspondence is seen 
only with the shift of three decimation levels between the 
two systems. Moreover, the critical value of the phase <j> 
setting the main peak at the origin is different for the 
models. In the QIM it is 1/4 (the same as at the lower 
band edge) as for the critical model the critical phase is 
3/8 which is the arithmetic mean of the critical phases 
at the upper and lower band edges. 

The Harper equation, which describes Bloch electrons 
on a square lattice, also models the isotropic XY model in 
a quasiperiodic field. H It is rather interesting that the 
problem for the Bloch electrons on a triangular lattice 
is related to the fully anisotropic Ising limit of the XY 
model. Although the correspondence between the TBM 
for the square lattice and the isotropic XY model is ex- 
act, relationship between the triangular lattice and the 
Ising model is through the asymptotic universal scaling 
properties of the models. 

Another interesting aspect of our studies is that unlike 
the periodic case, details of the lattice is relevant in de- 
termining the phase diagram and universality classes in 
incommensurate systems. For the square lattice, there is 
a phase transition from E to L phase while in the case of a 
triangular lattice, the C-L transition is observed instead. 
Therefore, for the square lattice, the weak coupling fixed 
point is trivial while for the triangular case, the weak 
coupling is described by the nontrivial fixed point and 
scaling ratios. Furthermore, even in the L phase, the 
scaling properties of the self-similar fluctuations in the 
square and the triangular lattices arc different. 

The crossover from E-L to C-L phase diagram, as we 
change from the square to the triangular lattice, can 
be understood due to the relationship between the 2- 
dimensional Bloch and spin problem. Previous studies 
Pj have shown that in the presence of anisotropy, the 
E and the L phases of the isotropic XY model are me- 
diated by a critical phase with self-similar characteris- 
tics. In the Ising limit, the width of the E phase shrinks 
to zero resulting in the C-L phase diagram. We like to 
mention that although the Harper equation with next- 
nearest- neighbor interaction [H 12 provides a means to 
study the cross-over from square to triangular lattice, 
understanding the E-L to C-L behavior is lot more com- 
plicated. 

The results of the triangular model may be of rel- 



evance in experiments on quantum dots. JlS) ] Further- 
more, since the one-dimensional QIM also describes the 
two-dimensional classical Ising model with modulating 
exchange along one direction, it could be realized in ex- 
periments on magnetic superlattices. Our studies pro- 
vide interesting mappings between the thermodynamical 
properties such as the specific heat and the susceptibility 
between these two models. 



ACKNOWLEDGMENTS 

The research of IIS is supported by a grant from Na- 
tional Science Foundation DMR 093296. IIS would like 
to acknowledge the hospitality of National Institute of 
Standard and Technology where part of this work is done. 
JAK has been supported financially by the Academy of 
Finland under contract 4385 and by the Niilo Helandcr 
Foundation. JAK is also grateful for the hospitality of the 
Research Institute for Theoretical Physics in Helsinki. 



[1] For a review, see J. B. Sokoloff, Phys. Rep. 126, 189 
(1985). 

[2] P. Christe and M. Henkel, Introduction to conformal 
invariance and its applications to critical phenomena. 
Springer, Berlin 1993. 

[3] J. A. Ketoja and I. I. Satija, Phys. Lett. A 194, 64 (1994); 
Renormalization approach to quasiperiodic quantum spin 
chains, to appear in Physica A. 

[4] F. H. Claro and G. H. Wannier, Phys. Rev. B 19, 6068 
(1979); D. J. Thouless, Phys. Rev. B 28, 4272 (1983). 

[5] D. R. Hofstadter, Phys. Rev. B 14, 2239 (1976). 

[6] J. A. Ketoja, Phys. Rev. Lett. 69, 2180 (1992). 

[7] J. A. Ketoja and I. I. Satija, Self-similarity and localiza- 
tion, (preprint) 1995. 

[8] E . Lieb, T. Schultz, and D. Mattis, Ann. Phys. (N.Y.) 
16, 407 (1961). 

[9] I. I. Satija, Phys. Rev. B 48, 3511 (1993); Phys. Rev. B 
49, 3391 (1994); I. I. Satija and J. C. Chaves, Phys. Rev. 
B 49, 13239 (1994). 

[10] We thank J.C. Chaves for confirming this. 

[11] J. A. Ketoja, I. I. Satija, and J. C. Chaves, Decimation 
studies of Bloch electrons in a magnetic field: higher or- 
der limit cycles underlying the phase diagram, Phys. Rev. 
B, in press. 

[12] J. H. Han, D. J. Thouless, H. Hiramoto, and M. 

Kohmoto, Phys. Rev. B 50, 11365 (1994). 
[13] R. R. Gerhardts, D. Weiss, and K. Von Klitzing, Phys. 

Rev. Lett. 62, 1173 (1989). 



3 



FIG. 1. The wave function in the subcritical (a), criti- 
cal (b), and in the supercritical (c) triangular model for the 
golden mean incommensurability at the upper band egde with 
= 0. Note that in the supercritical case the figure descibes 
the fluctuations \r)i\ and not \tpi\. 
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